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Abstract 


The occurrence of ductile fracture is the major limiting factor that determines the worka- 
bility of cold bulk metal form ing processes. Workability may be defined as the extent to 
which a material can be deformed without failure. Prediction of ductile fracture initiation 
will allow prior modification of the process which results in defect free final product with 
financial savings. 

Two common defects which occur in rolling are central burst and split-ends. The 
consequence of these defects are not only limited to less productivity and also causes 
damage to the rolls and mill accessories. Recent literature [Zhu and Avitzur, 1988; Avitzur 
et al., 1988; Turckzyn and Malinowski, 1996] as well as the current practice in the rolling 
industry shows that the successful avoidance of defect is largely a matter of empirical 
practice or based on the upper bound based fracture criterion. 

It is well known from literature [Dhar et al., 1996] that the ductile fracture occurs due to 
micro-void nucleation, growth and finally coalescence into micro-crack. Some experimental 
studies on extrusion and drawing revealed that the tensile triaxiality is the main factor 
contributing to the crack initiation. One can vizualize that the nature and magnitude 
of triaxiality generated in the deformation zone of plane-strain rolling process will be 
dependent on the various process variables. If the process variables are such that the 
tensile triaxiality arises in the deformation zone, then new micro- voids gets nucleated, 
existing voids grow and finally coalescence into micro-crack. Fracture occurrence can be 
prevented by choosing the proper process variables using mathematical modeling. Finite 
element method is the best modeling technique as it provides more accurate description 
of the deformation and stresses than other techniques. Hence, in the present work ductile 
fracture initiation is predicted using hydrostatic stress criterion [Reddy et al., 1996] and 
critical damage criteria [Dhar et al., 1996] along with the finite element results. 

Force, torque and pressure predictions of the present model are in good agreement with 
the available experimental [Al-Salehi et al., 1973 and Shida and Awazuhara, 1973] and 
analytical [ Dixit, 1997] results. Predictions of hydrostatic stress criterion for 
split ends are compared with the available experimental [Turckzyn and Malinowski, 1996] 



results and they are in good agreement. Parametric study is carried out to study the 
influence of process variables on the boundary between safe and unsafe zones for defects. 
It is observed that with increase in roll radius or decrease in friction the safe zone in 
increasing. 
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Chapter 1 

Introduction and Literature Review 


1.1 Introduction 

Rolling is one of the most important process in metal working industry, and at the same 
time it is one of the oldest metal forming processes. According to some estimations, 90% 
[ASM handbook, vol. 14] of the steel stock is processed at least once in a rolling mill. Plane 
strain rolling is the process in which the sheet metal is squeezed through the gap between 
the two counter rotating work rolls to reduce the thickness as shown in the figure 1.1. The 
rolling process is one amongst the most often used, economical and most efficient techniques 
available, to produce sheets, plates and strips. Rolling process because of its versatility, 
efficiency and its high production rate playing an important role in the manufacturing 
industry. Primary objectives of the rolling process are to reduce the cross section of the 
incoming material and at the same time improve its properties and to obtain the desired 
shape. 

Two common defects which occur in rolling are central burst (Fig. 1.2) and split ends 
(Fig. 1.3). The consequence of these defects are not only limited to less productivity 
and also causes damage to the rolls and mill accessories. In rolling process, without the 
complete knowledge of the influences of the variables such as friction conditions, material^ 
properties and workpiece geometry on the process mechanics, it will not be possible to 
predict and prevent the occurrence of defects, or to design the equipment adequately. 
It is almost impossible to find the exact analytical solutions of the equations governing 
real phenomenon by considering all aspects of process parameters and variables. Hence 
numerical simulation has become a major tool in modem metal forming research. Ability 
to predict the defect occurrence would allow prior modification of the process variables 
which eliminates the risk of failure occurrence, with obvious financial savings. 
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A number of approximate methods of analysis have been developed and applied to 
various metal forming processes. The most ■well known methods are the slab method, the 
slip line field method, upper bound method and more recently, the finite element method. 

Critical review of literature [Avitzur et al., 1988; Turckzyn and Malinowski, 1996] 
on defects in rolling process shows that the successful avoidance of defect is largely a 
matter of empirical practice or based on the upper bound based fracture criterion which 
is not accurate as described in a later section. It is well known from the literature [Dhar 
et al., 1996] that the ductile fracture occurs due to micro-void nucleation, growth and 
finally coalescence into micro-crack. Many empirical criteria have been proposed to predict 
the occurrence of ductile fracture in various metal forming processes. From literature, it 
is evident that the empirical criteria have had only limited success [Clift et al. 1990]. 
Investigations on ductile fracture reveal that it occurs due to micro- void nucleation, growth 
and finally coalescence into micro-crack. Quite a few void growth models have been applied 
to predict ductile fracture in certain forming processes. The porous plasticity models, which 
incorporates both void nucleation as well as void growth, have also been used in some cases. 
The main draw back of these models is that they do not treat micro-crack initiation as a 
void coalescence phenomenon. There are no attempts to predict the fracture initiation in 
rolling process by using the more realistic fracture criteria developed on the basis of void 
nucleation, growth and coalescence [Dhar et al., 1996, and Thomason, 1990]. Hence, it is 
necessary to predict and prevent the fracture initiation in rolling process by choosing the 
proper process variables using more realistic fracture criteria. 

The main aim of the present work is to predict and prevent the occurrence of central 
burst and alligatoring defects in plane strain rolling. Process modeling is performed using 
the finite element technique as it determines the process mechanics accurately compared 
to the other methods such as upper bound and slip-line methods. Two fracture criterion 
[Dhar et al, 1996 and Reddy et al., 1996] axe used along with the finite element results to 
predict the fracture occurrence in plane strain rolling process. Comprehensive parametric 
study is carried out to study the influence of process variables on the process parameters 
and defect occurrence. 
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1.2 Literature Review 

1.2.1 Deformation Analysis 

Rolling process has been analyzed extensively by the slab method, slip-line field method, 
upper bound method and the finite element method. 

1.2. 1.1 Slab method 

The slab method is known to be developed by Sachs and Siebel. In this method, the 
deformation zone is divided into a finite number of ’slab’ elements and an equilibrium 
equation for a typical slab is solved, subject to the given stress boundary conditions. The 
distribution of contact pressure at the tool-work interfax;e can be obtained based on a 
simplification of stress state. The assumptions on which this method is developed are: 

• The stresses on a slab surface, perpendicular to the flow direction, are working in 
principal directions and they are not permitted to vary on this surface. 

• Frictional stresses do not effect the internal stress distribution. 

Many researchers have been developed mathematical models to study rolling problem, 
using slab method. Starting from the pioneering work of von Karman [1925], Orowan 
[1943], Bland and Ford [1948], Bland and Sims [1953], Ford and Alexander [1963], Tselikov 
[1939], Stone [1953], Robert [1967] and Lenard and Roy Choudhuri [1984] and recently by 
Surya Narayana [1998]. A brief comparison of these models is given below. 

Orowan [1943], Bland and Ford [1948], Sims [1954] considered the roll contour to remain 
circular after the deformation and that the radius of curvature becomes larger than the 
original roll radius. They used Hitchcock’s [1935] method to calculate the deformed roll 
radius. Roy Choudhuri and Lenard [1984] pointed out that Hitchcock’s formula does 
not predict the deformed roll shape accurately and proposed an approximate method to, 
calculate the same. 

Pietrzyk and Lenard [1991] compared the predictive capabilities of the six [Orowan, 
1943, Sims, 1954, Bland and Ford, 1948, Ford and Alexander, 1963, Tselikov, 1939, Roy 
Choudhuri, 1984] models and along with experimental data. They observed that the 
Orowan’s [1943] technique proved to be reasonable in predicting roll forces and roll torques. 
The Bland and Ford [1948] analysis gave results very close to that of Orowan [1943]. They 
also observed that the inclusion of the experimentally determined data concerning friction 
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at roll strip interface enhances the ability of the model of Roy Choudhuri and Lenard [1984] 
in calculation of the roll forces. 

Lin et al. [1996] compared three one dimensional models qualitatively [Robert’s, 1967, 
Stone, 1953, Bland and Ford, 1948]. They reported that the Stone’s [1953] model, using 
Hitchock’s [1935] deformed radius equation gives a good estimation of roll separating force 
for strip reductions larger than 3 percent and less than 5 percent. They observed that 
when compared with other models, it is not recommended for reductions less than 3 per- 
cent. They also observed that Robert’s [1967] model using Robert’s [1967] deformed radius 
equation, can be used for predicting the roll separating force while rolling very thin strip 
with thickness less than 0.5 mm and reductions less than 3 percent. 

From the above discussion, it can be said that Robert’s [1967] method is efficient when 
compared with other methods as it considers strain rate effects, which other methods does 
not consider. The reason for this is the presence of the gradient of flow stress in the 
equilibrium equations in the other one dimensional methods mentioned above. Recently 
Surya Narayana [1998] developed a one-dimensional mathematical model, in which he 
eliminated the gradient of flow stress from the governing equations to incorporate the 
realistic material properties. The influence of strain rate in prediction of roll forces and 
roll torques was reported by Tseng et al. [1998]. 

1.2. 1.2 Slip-line method 

Slip line method is used to obtain the stress field in a forming problem by constructing 
slip lines. Originally, the method was proposed for strain analysis of rigid perfectly plas- 
tic materials. However, lately there have been some attempts to expand the method to 
axisymmetric problem as well as to incorporate strain hardening.’ 

Some researchers have analyzed the rolling process by slip-line method. Firbank and 
Lancaster [1965] proposed a slip line field for cold rolling with slipping friction. Their 
solution is for small rolls and a pass reduction of 20 percent. In slip line solutions peak 
pressures occur at the entry and exit, and there is considerable pressure drop some where in 
between those two points. New slip line solution for compression and rolling with slipping 
friction was presented by Collins [1969]. 

1.2. 1.3 Upper boimd method 

The upper bound method is a technique for determining an over estimation of the exact 
load, obtained from an assumed kinematically admissible velocity field which satisfies the 
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conditions of incompressibility and velocity boundary conditions. An extensive literature 
exists on upper bound analysis of metal forming processes and have been reviewed by 
many researchers [Avitzur, 1980, 1983; Bramley and Osman, 1992; Hill, 1950; Johnson 
and Mellor, 1973; Kudo, 1985; and Thomsen et al., 1965] and particularly on rolling by 
Avitzur [1988]. The upper bound method requires the guessing of admissible velocity fields, 
among which best one is chosen by minimizing the total energy. Information leading to a 
good selection of velocity fields come from experimental 'evidence and experience. 

A'vitzur’s [1964] upper bound solutions gives an expression for roll torque and the 
position of the neutral point, but does not provide roU force and pressure distribution 
which axe important design parameters. The major difficulty in using upper bound method 
is the selection of kinematically admissible velocity field since the accuracy of the solution 
depends on how close the assumed velocity field is, to the actual one. For the cases in 
which roll fiattening is significant and may affect torque, one has to compute roll force by 
some other method and calculate deformed radius of roll which is then to be substituted 
in the expression of the torque developed by A'vitzur [1964]. 

A'vitzur and Pachla [1986] developed a upper bound approach to plane strain deforma- 
tion of a rigid, perfectly plastic material in rolling process. They divided the deformation 
region into a finite number of rigid triangular bodies that slide -with respect to one another. 
Neighboring rigid body zones are analyzed in specific cases where the zones are (1) both 
in rotational motion, (2) one in linear and another in rotational motion, and (3) both in 
linear motion. Specific equations are presented that describes surface of velocity disconti- 
nuity (shear boundaries) between the mo-ving bodies, and the velocity discontinuities and 
shear power losses for each of the three cases. A'vitzur and Pachla [1986a] applied these 
equations drawing, extrusion and rolling. And exphcit equations are developed to calculate 
the surfaces of velocity discontinuity (shear boundaries), velocity discontinuities, and the 
upper bound on power for these processes. 

Avitzur et al. [1987] analyzed the strip rolling using the upper bound technique. The 
general solution given in Avitzur and Pachla [1986] is applied here to the problem of strip 
rolling. Two triangular velocity fields, one 'with triangles in linear rigid body motion and 
the other with triangles in rotational rigid body motion are developed. The total power 
is determined as a function of the four independent process parameters namely relative 
thickness, reduction, friction and net front-back tension. 
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1.2. 1.4 Finite Element Analysis 

Slab, slip-line and upper bound methods are useful in making the approximate prediction- 
s for the forming loads, the overall geometry changes in the deforming work pieces, the 
qualitative modes of metal flow and optimal process conditions. However, accurate de- 
termination of the effect of various process parameters on the metal flow became possible 
only after the flnite element method was introduced into the analysis of the metal forming 
processes. In flnite element method, the domain of interest is divided into a number of 
suitable elements and the governing differential equations of the system are reduced to 
algebraic equations by using appropriate approximations for the fleld variables over the 
elements. Almost all metal forming processes have been analyzed by the flnite element 
method due to its ability to handle a wide class of boundary value problems without re- 
striction on the work piece geometry. This can be achieved by selecting a proper mesh. 
Several review papers and text books have appeared outlining the state of art of the flnite 
element analysis of metal forming processes [Altan and Knoerr, 1992; Cser et al., 1993; 
Hartley et al., 1992; Herrmann, 1994; Altan and Vazquez, 1997]. 

Rigid-Plastic Finite Element Analysis 

The use of elasto-plastic formulation involves certain numerical difficulties. Since elastic 
effects are usually small, it is desirable to idealize the material as rigid-plastic, if it helps in 
avoiding the numerical difficulties. Zienkiwicz et al. [1978] considered the rolled material 
to be rigid-visco-plastic and incompressible. They simulated the friction at the roll-strip 
interface by introducing a thin layer of elements whose yield strength is assumed to depend 
on the coefficient of friction and mean stress. Mori et al. [1982] assumes that the rolled 
strip is made of a rigid-plastic, slightly compressible material. Using a constant frictional 
coefficient and rigid rolls they reported good agreement between the predicted and ob- 
served front end shapes of rolled aluminium strips. In their formulation, neutral point was 
determined by minimizing a certain functional. Li and Kobayashi [1982; 1984] also assumes, 
the existence of rigid-plastic materials and rigid rolls but they modeled the neutral point 
using the velocity dependent frictional stress. The authors included a comparison of their 
predictions with the results of Al-Salehi et al. [1973], and Shida and Awazuhara [1973]. In 
acknowledging that some discrepancies do exist, they attributed them to their use of rigid 
roUs and some imcertainty in modeling interfacial friction. 

Prakash et al. [1995] presented a FEM formulation in which the neutral point is found 
iteratively from the condition that the interfacial shear stress changes its sign at the neutral 
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point. This formulation does not incorporate roll deformation. Recently, Dixit and Dixit 
[1996] presented a rigid-plastic finite element formulation, to calculate the residual stresses. 
They applied fuzzy set theory to calculate material properties and friction coefficient. 
They applied this fuzzy set theory to avoid uncertain behavior in these parameters. They 
included anysotrophy in formulation, and also considered roll deformation. The front and 
back tensions also included in this formulation. 

Elasto-Plastic Finite Element Analysis 

Eulerian elasto-viscoplastic formulations have been reported by many authors. Dawson and 
Thompson [1978] have developed ’initial stress rate’ method. In this approach, the stress 
and velocity fields associated with the convective terms are assumed known and placed on 
the right-hand side of the finite element equations for momentum. An iterative procedure 
is then used during which these terms are corrected until there is no discrepancy between 
the new stress distribution and the precious one. The stress values are corrected until there 
is no discrepancy between the two successive iterations. Shimazaki and Thompson [1981] 
used the mixed formulation to solve for velocity, pressure and stress simultaneously. In 
their method, they decoupled the momentum and continuity equations from the constitu- 
tive equation. The former set was used to obtain velocity and pressure with stress assumed 
constant, the latter set was used to obtain stress with velocity and pressure held constant. 
Iterations between these two sets were done until convergence. Thompson and Berman 
[1984] employed the procedure developed by Shimazaki and Thompson [1981] for the anal- 
ysis of plane-strain rolling. However, rate of convergence was slow and results presented 
are without getting full convergence. EflFective stress contours are plotted, but individual 
stress components have not been presented. Thompson and Yu [1990] have presented a 
flow formulation to solve the rate equilibrium equation, in their steady state form, for ve- 
locity and stress. The authors have used the formulation for the analysis of strip drawing 
for a perfectly plastic material. It is noted that the rate of convergence deteriorates as the 
inelastic response becomes large. Effective stress contours are plotted without getting full 
convergence. No result is presented for residual stress distribution. 

Maniatty et al. [1991] have proposed an eulerian elasto-visco-plastic formulation in 
which the kinematics consists of a multipMcative decomposition of the deformation gradient 
into elastic and plastic components. They used Hart’s constitutive model, which contains 
an internal scalar variable, for analyzing plane strain rolling and extrusion processes. The 
system of equations is highly ill-conditioned and the algorithm lacks stability as the material 
becomes less rate sensitive. Maniatty [1994] modified the above formulation by using the 
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discontinuous pressure field by eliminating the pressure as a primary variable. He adjusted 
the mesh to align with the streamlines so that the constitutive equations can be integrated 
from node to node along the mesh lines and analyzed the rolling process to show the 
effectiveness of the algorithm. 

Malinowski and Lenard [1993] presented an elasto-plastic finite element formulation 
using the Prandtl-Reuss constitutive relations which was applicable for steady-state plane 
strain problems. In this formulation, the velocity field is first determined by minimizing 
the total power. The stresses are then obtained by integrating the Jaumann stress rate 
along the pathlines and minimizing a functional which relates the Jaumann rate to the 
material derivative of the stress using the finite element method.' Pathline integration 
gives unrealistic results due to non-satisfaction of the equilibrium condition. In the finite 
element procedure of finding stresses, they used two types of shape functions (parabolic and 
Hermitian) and the results differ significantly from each other. The formulation assumed 
the poisson’s ratio to be a continuous material function defined for both elastic and elasto- 
plastic behaviors. Malinowski [1993] used the methodology developed by Malinowski and 
Lenard [1993] for determining residual stress distribution in plane strain cold rolling. 

Recently, Dixit and Dixit [1997] presented elasto-plastic formulation for rolling problem. 
The formulation presented here is similar (but not exactly the same) to that of Shimazaki 
and Thompson [1981], Thompson and Berman [1984], the main difference being in the 
constitutive equations. In Shimazaki and Thompson [1981], Thompson and Berman [1984], 
they considered whole domain as elasto-plastic and the same constitutive equation is used 
throughout the domain, where as in Dbdt and Dbdt [1997] different constitutive relations 
are used for the elastic and elasto-plastic regions. 

1.2.2 Ductile Fracture Criteria 

Ductile fracture is often a limiting factor in many metal forming processes. It is well known 
that ductile fracture occurs due to micro- void nucleation, growth and finally coalescence, 
into micro crack. Many ductile fracture criteria have been proposed and applied to predict 
ductile fracture in plastically deforming metals including metal forming processes. This 
section reviews some of the published ductile fracture criteria. 

The published ductile fracture criteria can be broadly classified into two groups. They 
are 

• Models based on microscopic observations about void nucleation, growth and coales- 
cence 
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• Empirical and semi-empirical models 

Three broad approaches have emerged which try to predict ductile fracture (micro-crack 
initiation) on the basis of void nucleation, growth and coalescence. They are 

• Models based on porous plasticity 

• Models of void nucleation, growth and coalescence 

• Continuum damage mechanics models 

1.2.2. 1 Models of Porous Plasticity 

In these models, the material with voids is idealized as porous material. Thus, its con- 
stitutive equation is derived form the plastic potential of porous materials. In a porous 
plastic material, there is a possibility of dilatation because of which the yield surface, does 
not remain an infinitely long cylinder like that of an incompressible material but is capped 
by elliptical surfaces. 

Gurson's Model 

Based on the Berg’s [1970] model of dilatational plasticity, Gurson [1977] proposed a plastic 
potential function for a porous solid with randomly distributed voids of volume fraction 
Vf in the form 

$ = ^ + 2V>cosli(^)-(l + V/)=0 (1.1) 

where ap and are the generalized and hydrostatic stresses of the porous aggregate and 
a is the generahzed stress of the void free (incompressible) matrix. 

In this model, the growth rate of void volume fraction is considered partly due to the 
growd;h of the existing voids and partly due to the nucleation of new voids, i.e., 

Vf — Vf— growth "k Vf— nucleation (l'2) 

The rate of change of the void volume fraction due to the growth of existing voids is related’ 
to the mean strain rate e,-,- by the relation 

Vf.growth = (1 - Vf) ea (1.3) 

In general, void nucleation rate depends on both the generalized strain rate as well as 
the rate of increase of hydrostatic stress. In this model, void nucleation rate is assumed to 
depend only on the generalized strain rate by the relation 

Vf -nucleation ~ (I-^) 
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where A is a constant. 

Gurson’s criterion [1977] was modified by Tvergaard [1981] and Needleman and Tver- 
gaard [1984] by introducing some constant (51) to bring the predictions of the model into 
closer agreement with numerical analysis of a periodic array of voids. The modified crite- 
rion is given by 

4 = ^ + 2ftV;Ccsh (5^) - (1 + (51^/)=^) = 0 (1.5) 

where qi is a constant. 

In this model, ductile fracture is regarded as the result of an instability in the dilata- 
tional plastic flow field localized in a band (called the shear band). The fracture criterion is 
represented as a graph of critical localization strain versus the critical void volume fraction 
with strain hardening exponent as a parameter. 

Oyane’s Model 

Oyane [1972] and Oyane et al. [1980] used the plasticity theory for porous materials 
to propose a ductile fracture criterion indicating that micro crack initiates whenever the 
volumetric strain reaches a material dependent critical value. For pore free materials, 
assuming that the material obeys porous plasticity equations after the initiation of voids, 
Oyane obtained the criterion for ductile firacture as 

( 1 . 6 ) 

where fp is a function of relative density Pr which is defined as the ratio of the apparent 
density of the porous material to the density of the pore free matrix, is the volumetric 
strain at fracture and n and C are constants. The constant C is given by 

where Cj is the generalized strain at which voids get initiated, e/ is the generalized strain' 
at which fracture occurs and A is a material constant. By assuming that e, = 0, the above 
criterion reduces to a simple form 

( 1 . 8 ) 

For porous materials, they modified the criterion by including the relative density term 
into the integral. Thus,' 
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Here is the generalized flow stress of the porous material, B is a constant and po is 
the initial relative density of the porous material. A method for estimating values of the 
material constants using the compression test has been provided by them. 

The major limitation of the above two criterion is that they do not model ductile 
fracture (i.e., micro-crack initiation) as a void coalescence phenomenon. 

1.2. 2. 2 Models of Void Nucleation, Growth and Coalescence 

The void nucleation model of Goods and Brown [1979], the growth models of McClintock 
[1968] and Rice and Tracey [1969] and void coalescence model of Thomason [1990] along 
with some additional models on void nucleation and growth are discussed in this section. 

Models on Void Nucleation 

Goods and Brown [1979] developed a stress based void nucleation model based on condition 
that the voids nucleate by decohesion of particles whenever the normal stress reaches a 
critical value (cTc) at the particle/matrix interface. They obtained a relation between the 
void nucleation strain (e” ) and hydrostatic stress (a^) as 

e" = Arrp((7c - Gmf (1.10) 

where K is a material constant dependent on the particle volume fraction and Vp is the 
particle radius. Experiments on Fe — Fe^C system [Goods and Brown, 1979] confirm the 
hnear relationship between the nucleation strain and particle radius. This relation is for 
small spheroidal particles of radius less than or equal to lyum. 

For large size particles, the approximate analysis of Argon and Im [1975] gives the 
condition for the void nucleation by decohesion as 

^ + crm = CFc ( 1 - 11 ) 

For a given problem, this condition can be expressed in terms of the nucleation strain. 
Note that this relation is independent of particle radius. 

Gurland [1972] proposed a nucleation model based on the experiments conducted by 
him on 1.05%C spheroidal steel. Experimental observations reveal that the void nucleation 
occur by cracking of cementite particles. Voids nucleate at all strain levels depending on 
the size, shape and orientation with the maximum principal stress direction. It is observed 
that the volume fraction of broken particles is a linear function of the plastic strain. 
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Void Growth Models 


Once the nucleation of micro-void takes place either by decohesion or cracking of a sec- 
ondary particle or inclusion, the resulting stress free surface of the void produces a localized 
stress and strain concentration in the plastically surrounding matrix. With continuing plas- 
tic flow of the matrix, the void undergoes volumetric growth and change of shape. Some 
of the published void growth models are discussed here. 

McClintock [1968] proposed a void growth model for two-dimensional plane strain prob- 
lems considering a single elliptic void. He assumed that the major and minor axes of the 
void coincide with the principal stress directions. He obtained the following closed for- 
m expression for void growth in a rigid work-hardening material whose generalized stress 
relation is a ^ K^\ 



2(1 — n) 2 (7 


+ 


+ (-b 
2 


( 1 . 12 ) 


where and R^h E^re the current and initial mean radii of the hole respectively, Uq, and aj 
are the principal stress and ta and ej are the principal strain components along the major 
and minor axes respectively. Fracture was assumed to occur at the point where a growing 
void touches the cell boundary. Dung [1992a, 1992b] modified and extended McClintock’s 
model [1968] to analyze the growth of circular/ellipsoidal void in three-dimensional prob- 
lems. 


Rice and Tracey [1969] considered a single spherical void of initial radius Rov hi a remote 
uniform strain rate field and remote stress field Cij in a plastic material. They obtained 
the following expression for the rate of change of the radii (i?A;v) of the void in the principal 
strain rate directions: 


Rkv — {(1 + -FjCjc -f €fr}i?rn ean (1.13) 


where 
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F 


H 




3 for linear hardening and for low values of cr^ for non hardening 
1 for high values of cr^ for non hardening 
0.75^ for linear hardening 

0.558 sinh + 0.008i/cosh for non hardening, 

generalized strain rate 
principal strain rates 
~J^^> Lode variable. 


Note that an initial spherical void grows into an ellipsoidal void of principal radii 
-^Sv Thomason [1990] integrated the above equation by assuming that the 
principal axes of strain rates remain fixed in direction throughout the strain path and 
obtained the following expressions for the principal radii of the void; 


where 


— A -\- 


i?2v = 

Rsv = 


i"- 

A + 


35(1 + 

2 + 3 


(1.14) 

Bv \ 

Rov 

(1.15) 


B{v - 3 ) ^ 

2Vj^2 -+3 ; 

1 Rov 

(1.16) 




B = 


3 + 1 ^ 
1 + F 
H 


(. 4 - 1 ) 


and e( is the integral of the largest principal strain rate. Thomason [1990] used the above 
expressions in the derivation of the void coalescence condition. 


For an array of void nucleating particles of diameter Dp and spacing ds, setting the 
initial void radius as Dp/2, eqns. (1.12) and (1.13) leads to the following expression for 
fracture strain: 


In {da /Dp) (1 -n) 

sinh [(1 - n)((7a + ot/) / {sd / \/Z)\ 


[McClintock, 1968] 


(1.17) 


_ _ In {da /Dp) 

0.28 exp (1.5(jjfcfc/CT) 


[Rice and Tracey, 1969] (1-18) 
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Some research workers have used the above expressions for the prediction of ductile fracture 
with the help of experimentally determined parameters. The major limitations of the 
approach are that it ignores the effects of void nucleation (as the expressions are derived 
from a pre-existing finite size void) and void coalescence. Thus, ey is usually overestimated. 

Void Coalescence Model 

Thomason [1990] modeled void coalescence phenomenon as plastic instabihty du'e to neck- 
ing of the inter void matrix. The sufficient condition for plastic instability of the inter- void 
matrix, as given by him, is 

OnAn - (Ti = 0 (1.19) 

where An is the area fraction of the inter-void matrix perpendicular to the direction of 
maximum principal stress ax and an is the plastic constraint stress. He considered a 
geometrically equivalent square prismatic void and used the upper bound method to obtain 
an expression for the plastic constraint stress in terms of the void dimensions, inter- void 
spacing and the yield stress. 

Thomason's Fracture Criterion 

Thomason [1990] combined the results of Goods and Brown [1979] on void nucleation. 
Rice and Tracey [1969] on void growth and his own on void coalescence to arrive at a 
fracture criterion. He used expressions (eqn. 14 to 16) for the void dimensions to express 
the void coalescence condition (eqn. 19) in terms of the void growth strain e? and the 
hydrostatic stress am- By superposing this condition on the void nucleation relation (eqn. 
10), he obtained the fracture criterion as a graph of fracture strain (e/ = €" -H ef) versus 
the triaxiality (i.e., the ratio of hydrostatic stress to the generalized stress). 

The main drawback of this approach is the use of expressions (eqn. 14 to 16) for void 
growth. There expressions are based on integration procedure which assumes that the 
principal directions of strain rate remain fixed throughout the strain path. This is, in 
general, true only of the case of small strain and rotation. 

1.2. 2. 3. Continuum Damage Mechanics Models 

In these models, material behavior is represented by a plastic potential which includes 
damage as an internal variable. The damage variable quantifies the intensity of voids, which 
can be identified as either the void volume fraction or the area void fraction. The basis of 
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continuum damage mechanics models rests on the theory of continuum thermodynamics. 
In this section some of the continuum damage mechanics models are discussed. 

Lemaitre [1985] proposed a continuum damage mechanics model for void growth in 
elasto-plastic material with area void fraction as the damage variable using concepts of 
effective stress and strain equivalence. His model is based on a simple damage growth law 
in which the damage rate is linearly dependent on the plastic strain and the thermodynamic 
force corresponding to the damage. It does not account for void nucleation. Note that, as 
far as void coalescence is concerned, it has to be incorporated as an additional condition in 
terms of the continuum parameters. This condition, which serves as a fracture criterion, 
has to be based on an appropriate micro model. Thomason’s limit load model seems to be 
a good candidate for this purpose. 

More recently Dhax [1995] and Dhar et al. [1996] extended Lemaitre’s [1985] continuum 
damage mechanics model to incorporate a damage growth law which accounts for both the 
void nucleation as well as void growth. The law is based on the experimental results of 
Le Roy et al. [1981]. They also modified Thomason’s [1990] void coalescence condition to 
make it suitable for the case of large strain and/or rotation. They combined the extended 
continuum damage mechanics model and modified void coalescence condition to arrive at a 
criterion for prediction of micro crack initiation. The salient features of the above approach 
are discussed below. 


In this model, damage D is identified as an area void fraction at a point in a plane. 
Thus, it is defined by 

D = ^ (1-20) 

where AA is the infinitesimal area around the point in some plane and AA^ is the area of 
the voids in the plane containing AA. The conjugate variable corresponding to D is -Y, 
the rate at which the elastic energy release during the damage growth at constant stress. 
For an isotropic material, -Y is given by [Lemaitre, 1985] 


-y = 


2E{1 - Df 



( 1 . 21 ) 


where 


/(^) = + 3(1 _ 


Here, E is the Young’s modulus, v is the Poisson’s ratio and 
the damage growth law is given by 


( 1 . 22 ) 

is called the triaxiality. 


t) = (fli + <X2D){—Y)(. 


(1.23) 
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where oi , CI 2 Cd &re material constants. Here the first term, which is independent of 
-Y, represents the damage evolution due to void nucleation while the second term, which 
is dependent on -Y, represent the evolution of damage due to void growth. For SAE 
1090 steel, Dhar et al. [1995] determined the constants fitting eqn. (1.23) through the 
experimental results of Le Roy et al. [1981]. They obtained the following values; 


Cd = 1.898X10-2 

(1.24) 

oi = 9,8X10"^ MPa-^ 

(1.25) 

02 = 1.84 MPa-^ 

(1.26) 


While getting these constants, they used Bridgman’s relation [1964] to express the triaxi- 
ality as a function of strain. 

Dhar et al. [1996] used finite strain expressions for the void dimensions and inter- void 
spacing to express Thomason’s void coalescence criterion (eqn. 1.19) in the following form: 


CTl 


0.1 + 


1.2 

{1 -ea:p(-e/2)}°-^ 


exp(—€) cr = 0 


(1.27) 


where e, the generalized strain, is the integral of t along the path lines. To find the 
critical value of damage parameter (Dc), they apphed this criterion to various geometries 
and loading conditions using the elasto-plastic finite element analysis. They observed that 
Dc is independent of geometry or loading and here that can be used as a material property 
for the prediction of micro-crack initiation. They reported that the value of Dc is 0.05 for 
SAE 1090 steel. 


1.2. 2.4. Empirical and Semi-Empirical Models 

In the absence of reliable quantitative models for incorporating the phenomena of void 
nucleation, growth and coalescence in materials undergoing large plastic deformation in 
metal forming processes, many empirical relations based on some phenomenological models 
have been proposed. This section reviews some of them. 

Frudenthal’s Model 

Preudenthal [1950] postulated that generalized plastic work per unit volume is the critical 
parameter and is expressed as 

= (1.28) 
Jo 
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Cockcroft - Latham's Model 


Cockcroft and Latham [1968] modified the generalized plastic work criterion to take care of 
the effects of change in the neck geometry which is observed in tensile tests. The criterion 
is given by 

Jo ''**■= ( 1 - 29 ) 

Here ai is the maximum normal stress. 


Oh’s Model 

Oh et al. [1979] modified Cockcroft and Latham criterion [1968] by considering the relation 
of the maximum tensile stress to the generalized stress. This criterion is expressed as 

= (1.30) 

Jq a 

The main drawback of the above three criteria is that the effect of hydrostatic stress is , 
not considered. 


Brozzo's Model 


Brozzo et al. [1972] proposed an empirical modification of Cockroft and Lathat’s model 
to consider the effect of hydrostatic stress explicitly. The criterion considers Hydrostatic 
stress is given by 

ri 9/t, 

de = CA (1.31) 


I 


0 3(cri 


0 


Norris's Model 

Norris et al. [1978] developed a fracture criterion based on their experimental results and 
finite difference results of several test geometries and expressed it as 

where C is a constant. 

In the above equations Ci , . . . Cs are the critical material dependent values accumulated 
over the strain path to fracture. 



1.2. 2. 5 Hydrostatic Stress Criterion 


Some experimental studies on extrusion revealed that the tensile triaxiality is the main 
factor contributing to the crack initiation. One can vizualize that the nature and mag- 
nitude of triaxiality generated in the deformation zone of plane-strain rolling process will 
dependent on the various process variables. If the process variables are such that the ten- 
sile triaxiality arises in the deformation zone, then new micro- voids gets nucleated, existing 
voids grow and finally coalescence into micro-crack. The work of Clift et al. [1990] shows 
that the value of hydrostatic stress at fracture initiation in extrusion is close to zero. Based 
on this observations, Reddy et al. [1996] proposed a simple criterion given as 
’whenever hydrostatic stress at a point on the centerline in the deformation zone becomes 
zero and is compressive elsewhere, there is fracture initiation leading to central burst’. 
One of the main advantage of this criterion is that it does not involves any critical material 
parameter which needs to be determined experimentally. 

1.2.3 Defects in Rolling Process 

DifiFerent type of defects are observed in industrial metal forming processes which have 
been tabulated by Johnson and his co-workers [Johnson and Mamalis, 1977; Mamalis and 
Johnson, 1987; Johnson, 1991]. Two common defects which will occur in rolling process 
are central burst and alligatoring (split-end) [Avitzur et al., 1988]. 

Criterion for prevention of central bursting and split ends in plane-strain rolling were 
derived by Avitzur and his coworkers [Zhu and Avitzur, 1988; Avitzur et al., 1988] by ex- 
tending an earlier upper bound solution [Avitzur, 1964] for strip rolling of a rigid perfectly- 
plastic Mises material. Their criterion for defect prediction is that for a given set of process 
parameters the power consumption for defect-free flow is higher that that for a defect flow. 
Turckzyn and Malinowski [1996] used the same criterion as above to predict the central 
burst and split ends in rolling. They also conducted experiments to predict occurrence of 
split end defects in rolling. 

It is well known that the tensile triaxiality is the main factor contributing to the crack 
initiation in metal forming processes. Note that, all the attempts made for predicting the 
central bursting and split ends in rolling [Zhu and Avitzur, 1988; Avitzur et al., 1988; Tur- 
ckzyn and Malinowski, 1996] are based on kinematically admissible velocity field (Upper 
bound analysis) but not on the basis of statically admissible velocity field. The circum- 
stances which may lead to deduction of a greater driving force for the process involving a 
defect must arise because the upper bound to the driving force with a defect is sufficiently 
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above the correct value to exceed the greater correct value corresponding to defect-free 
product. Use of lower bound theorem [Lee and McMeeking, 1978] indicates that the actual 
driving force for the cracked workpiece provides a lower bound for the uncracked case. 
Thus the correct driving force for the uncracked case is always greater than or equal to 
that for the cracked case. Hence, the inequality in the defect prediction criterion used by 
Avitzur and his coworkers [Zhu and Avitzur, 1988 and Avitzur et al., 1988] and Turckzyn 
and Malinowski [1996] implies that the defect always occurs. 

1.3 Scope and Objective of the Present Work 

Critical literature review reveals that there have been very few attempts to predict the 
occurrence of central burst and alligatoring in rolling. Recent literature [Avitzur et al., 1988 
and Turckzyn and Malinowski, 1996] as well as the current practice in the rolling industry 
shows that the successful avoidance of defect is largely a matter of empirical practice or 
based on the upper bound based fracture criterion. Present work aims at prediction and 
prevention of central bursting and split ends in plane-strain rolling by choosing the proper 
process variables through finite element based process modeling as it determines the process 
mechanics accurately compared to other methods such as upper bound method and slip 
line method. The various rolling process variables which will effect the defect occurrence 
are relative reduction in thickness, front and back tensions, friction, material properties, 
roll velocity and radius. 

The main objective of the present work is to predict and prevent the occurrence of 
central burst and alligatoring in rolling. The above objective is achieved through the 
following steps: 

• Development of Eulerian rigid-plastic finite element formulation with isotropic hard- 
ening and implementation 

• Validation of finite element predictions with available experimental and numerical 
results 

• Identification of proper ductile fracture criteria 

• Implementation of ductile fracture criteria in the computer program for fracture 
prediction 

• Validation of predictions regarding fracture occurrence 

• Study the effect of process variables on fracture occurrence 
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1.4 Organization of the Thesis 

The thesis chapters are organized in the following manner: 

Chapter 1 discusses the importance of prediction of defects in rolling process. Critical 
literature review on rolling related to deformation analysis by various methods and ductile 
franture is presented. Based on the literature review the scope and objective of the present 
works is presented. At the end, the way in which the thesis is organized is presented. 

Chapter 2 is concerned with the Galerkin finite element formulations to obtain the de- 
tailed deformation field in plane-strain cold rolhng. A rigid-plastic formulation is presented 
where the elastic effects are considered to be negligible. In general the generalized stress in 
the cold rolling material is a function of strain. Therefore, a finite element formulation for 
the determination of generalized strain distribution is presented. To predict the occurrence 
of damage by the critical damage criterion of Dhar et al. [1996], distribution of damage 
variable D over the domain is required. For this distribution of D, a Galerkin finite element 
formulation is presented which uses the damage law employed by Dhar et al. [1996]. To 
predict, the occurrence of fracture by empirical criterion of Oyane et al. [1980], a Galerkin 
finite element formulation is presented to determine the distribution of damage variable. 
At the end, procedure to obtain the secondary variables such as roll torque, roll force and 
interfacial normal pressure is presented. 

In chapter 3, first different fracture criteria used are presented in section 3.2. In section 
3.3, predictions of present model are compared with the available experimental and analyt- 
ical results. Parametric study is carried out to study the influence of process variables on 
hydrostatic stress distribution, damage distributions of critical damage criterion [Dhar et 
al., 1996] and Oyane’s criterion [Oyane et al., 1980] and the same is presented in sections 
3.4.1 to 3.4.3. The effect of process variables on the safe and unsafe zones for fracture is 
presented in section 3.4:4. 

In chapter 4, conclusions and scope for future work is presented. 
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Chapter 2 

Finite Element Formulation 

2.1 Introduction 

An important step in the prediction of defects is to develop a fairly exact mathematical 
model of the process for the determination of the deformation and stress fields. An accurate 
determination of these fields became possible after the finite element method was developed 
for the metal forming processes. 

In this chapter, Galerkin finite element formulations are presented to obtain the detailed 
deformation field in steady-state plane-strain rolling process. In section 2.2, a rigid plastic 
formulation is discussed where the elastic effects are considered to be negligible. In general, 
the generalized stress of the material is a function of strain for the case of cold rolling. 
Therefore, one has to know the strain distribution for accurate estimation of the generalized 
stress. The procedure for determination of strain distribution is presented in section 2.3. 
To predict the occurrence of firacture by the more realistic criteria, distribution of damage 
variable over the domain is required. Galerkin finite element formulations are presented in 
sections 2.4 and 2.5 to determine damage distribution. In section 2.4 damage formulation 
is for critical damage criterion [Dhar et al., 1996] and in sections 2.5 for Oyane’s criterion 
[Oyane et al., 1980]. The procedure for obtaining the secondary variables is presented at 
the end of this chapter. 

2.2 Rigid- Plastic Formulation 

In bulk metal forming processes such as rolling, extrusion, drawing etc., the plastic strains 
are very large compared to elastic ones. Hence elastic strains can be considered to be 
insignificant as far as the deformation phenomena are concerned. When the elastic effects 
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are neglected, then the material behavior is essentially that of a fluid i.e., material incapable 
of sustaining any deviatoric stress, without motion. Under the assumption of negligible 
elastic strains, the material can be idealized either as rigid-plastic or rigid-viscoplastic. 
Deformation of such materials can be conveniently analyzed by Elurian or flow formulation 
[Zienkiewicz, 1984] and it is with this formulation for plane-strain rolling that this section 
is concerned. 

2.2.1 Constitutive Equation of Rigid-Plastic Materied 

The system of equations describing the behavior of a solid under various types of response 
are referred to as the constitutive equations. Steady-state forming processes such as rolling, 
extrusion, drawing etc., can be conveniently analyzed by using an Eulerian formulation. 
For such a formulation, the measure of deformation is the strain rate tensor Cy, which is 
expressed as 

ey - 2 

where Vi is the velocity component and Xj represents the coordinate system. In order to 
express the constitutive equation in a convenient form, the stress tensor Uy is decomposed 
as 

(Tij = -p5ij + c/ij] p=-^ (2.2) 

where p is the (compressive) hydrostatic part, 5y is the Kronecker’s delta and is the 
deviatoric part. Henceforth, the hydrostatic part of Oy vdll be denoted by -p and not by 

The constitutive law for rigid-plastic/ viscoplastic material relating the deviatoric stress 
tensor and the strain rate tensor is expressed as 

(Vij - 2/i€y (2.3) 

For the material yielding according to von-Mises criterion, the Levy-Mises- coefficient jj. 


is given by 



IX- 

cy 

(2.4) 

where the generalized stress and strain rate are defined by 


and 

/a 

1 

(2.5) 


[2. . 

1 

(2.6) 


^ doCj doCi ^ 


( 2 . 1 ) 
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In general a depends on e, e and temperature T i.e., 

Here, the generalized strain e is defined as 


e = / kdt 


where the integration is to be carried along the particle path. 


2.2.2 Governing Equations 


In general, the material behavior of the deforming metal is governed by the equations of 
conservation of mass, conservation of momentum and conservation of energy. The present 
formulation deals with the steady-state plane-strain cold rolling where the forming speeds 
are low and therefore the temperature rise is small. Thus, the process is considered to 
be isothermal. Therefore, the equation of conservation of energy is not required for the 
analysis. 

In Eulerian reference frame, the momentum equations for a steady-state plane-strain 
rolling process (neglecting the body forces) have the following form: 


.. dv,. 


I ^^xy 

dx dy 


^ ^ dy dx dy dy 


( 2 . 10 ) 


where p is the density, Vx and Vy represent the velocity components with respect to the 
coordinates x, y (Fig. 2.1), p is the hydrostatic part and is the deviatoric stress part 
on the stress tensor Oij. In plastic deformation, volume remains constant. Therefore, the 
hydrostatic part of stress is determined from the continuity equation 


^xx d* — 0 


( 2 . 11 ) 


The relations between strain rate and velocity components for plane-strain problem are 
given by 


( 2 . 12 ) 


^yy — 


~ n 


dy 

2^ dy dx 


( 2 . 13 ) 

( 2 . 14 ) 
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2.2.3 Boundary Conditions 


In the present work, plane-strain strip rolling processes is considered for the analysis. The 
boundary conditions for the same are discussed in this section. 

The boundary conditions specified on any boundary can be classified into three cate- 
gories. They are velocity boundary conditions (specification of velocity components), stress 
boundary conditions (specification of traction components) and mixed boundary conditions 
i.e., combination of one velocity and one stress boundary condition. 

The domain along with co-ordinate system for rolling process is shown in Fig. 2.1. 


Entry and Exit boundary conditions (AB and EF) 


The entry and exit boundaries AB and EF respectively are chosen sufficiently away from 
the deformation zone so as to take the advantage of uniform conditions prevailing there. 
Therefore, the conditions on the boundaries AB and EF become 


Vx = Kntry, Fj, = 0 on AB (2.15) 

14 = "K, = 0 on EF (2.16) 

where V^ntry and Vtxu are respectively the inlet and exit velocity related by 


Gentry — ^exit (1 ^) 


(2.17) 


where r is the reduction ratio given by r = 

In a rolling problem inlet and exit velocities are not specified, they can be calculat- 
ed provided the position of neutral point and roll velocity using the principle of volume 
constancy. 

If front and back tensions are zero, surfaces AB and EF are stress free, then the following 
boundary conditions can also be applied. 


Fj, — 0; tx = 0 on AB and EF (2.18) 

Free Surfaces (BC and DE) 


The boundaries BC and DE are free surfaces, therefore, both t^ and ty are zero on BC 
and DE. It is also possible to use the condition Fj, = 0. In the present work, the following 
boundary conditions are used: 

Fy = 0 and tx = 0 ion AB and CD (2-19) 
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The Plauie of Symmetry (AF) 

The boundary AF is the plane of symmetry, therefore, the following conditions are used: 

tx = 0 and = 0 on AF (2.20) 

The Roll- Work Interface (CD) 

On the roll-work interface, the normal velocity component at any point must be zero. Thus 

= 0 => -f- 14 tan (;6 = 0 (2-21) 

where (f) is the angle made by the outward unit normal to the interface at that point with 
the y-axis. 

The second boundary condition is specified in terms of a friction condition. In the 
present model Coulomb’s friction model is used. When Coulomb’s law is used, normally a 
constraint is placed on the maximum value of the shear stress. Thus, 

\is\ = f\tn\ when |t,| < ^ (2.22) 

and 

Here and tn are the tangential amd normal stress components along the interface and / 
is the coefficient of friction assumed to be constant along the interface. 


2.2.4 Galerkin or Weak Formulation 


In the present work, mixed pressure-velocity formulation is used. The Galerkin weight- 
ed residual method is used to reduce the governing diflferential equations of the problem 
into algebraic equations by using appropriate weighted functions for the momentum and 
continuity equations. 

Let Vx, Vy and p be the functions which satisfy all the essential boundary conditions 
exactly. Then I 4 , Vy and p constitute a weak solution if the integral equation. 



+ ^yv)^p + 





dy dx 


dx 






dxdy = 0 


(2.24) 
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is satisfied. Here Wp, Wx and Wy are the approximate weight functions which satisfy the 
homogeneous versions of the boundary conditions and represents the area of the domain. 

Performing the integration by parts, on second and third parts of the above equation 
and using eqn.(2.3), the following weak form is obtained 


l^[h + 12] dx dy- dTi - ^ J4 dTj = 0 


(2.25) 


where 


h — {^xx + ^yy) 


(2.26) 


dVx , ,, dV, 


dV. 




-pfe.(M') + imiW)) + + 2e^€^(W')}} (2.27) 


h = txW'. 


(2.28) 


I 4 = tyWy 


(2.29) 


and Pi and r 2 are the boundaries where traction components tx and ty are specified. The 
terms ixx{W), iyy{W) and ixy^W) are defined as 


exx(W) = 


iyy(W) = 


dx 

dW, 

Sy 


, = 1(^ + ^) 
2^ dy dx^ 


(2.30) 

(2.31) 

(2.32) 


For finite element formulation, it is convenient to express integral (2.25) in the matrix form 
as 


L 


( 8^ ^ 


{mY{e}Wp + p{Vx{Wx Wy} 


dVr 

9y 


} + Wy}\ }} - p{Kw)Y M 

eVy 8 Vy 

< dx ' I 5v 


+2)u{e(lT)}^{€}] dxdy - txWxdT^ - tyWydr 2 = 0 (2.33) 
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where 


and 


/ • N 


r ^ 

\ 

^XX 


■dx 


% 

► = 






dx ) J 





dx 

= < 

iyylW) 

= < 

dWy 

dy f 


. V2i^{W) , 




r 1 

\ 

{m} = 

1 



1 0 

* 


( 2 . 34 ) 


( 2 . 35 ) 


(2.36) 


2.2.5 Finite Element Approximation 


In contrast to the traditional Galerkin weighted residual method where a global approx- 
imation is chosen for the field variable, in finite element method the domain of interest 
is discretized into small elements and the integrals are evaluated over each element after 
suitably selecting the approximation for the field variable over each element and then the 
elemental expressions are assembled appropriately to get the global expressions. 


From eqn. (2.25), one can see that the convergence and completeness criterion are 
satisfied if Vx and Vy are chosen to be bilinear and p is chosen to be constant. However, 
one can use higher order approximation like biquadratic approximation for Vx and Vy and 
bilinear approximation for p. In the present work, biquadratic approximation for velocities 
and bilinear approximation for pressure are used (Fig. 2.2). Thus, the approximations for 
Vx and Vy are 



' N, 

0 •• 

1 

O 

[Vyl~ 

0 

Ni-- 

• • 0 iVg J 


(^ 1)9 

(^y)9 




(2.37) 


In Galerkin formulation, the weight functions corresponding to velocities are approximated 
using the same shape functions as that for velocity. Therefore, 

{W,) = { ^' } = (2-38) 
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The vectors {V‘} and contain the nodal values of velocities and the corresponding 

weight functions, respectively. The approximation for p gives 


Pi 

P\ 


p = { JV| JVJ iVT } , , 


P\ 


= mif} 


(2.39) 


The weight functions corresponding to pressure are approximated using the same shape 
functions as that for the pressure. Therefore, 


W, = \N,]{W\] (2.40) 

The vectors {p*} and contains the nodal values of pressure and the corresponding 

weight functions, respectively. The expressions for the biquadratic {Ni) and bilinear {N^i) 
shape functions are given in any standard text book [Reddy, 1993] which deals with the 
finite element method. 

In the present work, the geometry is also approximated by the same shape functions as 
that used for velocities. Therefore, 


I = [JV]{a:'}; y = [iV]{y'} (2.41) 

The vectors {r®} and {j/®} contain the nodal values of x and y coordinates respectively 
and {A^} is the one-dimensional array of quadratic shape functions as 

{N} = {Ni N2 • • • iVs } (2.42) 


To evaluate the integrals over boundaries Fj, an approximation for the weight func- 
tions over these boundaries is needed. The approximation consistent with the interior 
approximations (eqn. 2.38) is 


W* = 


JV* iV* Nl] 


'''vx 

"vx 
W3 
'^^vx J 




and 


lV,= [iV{ JV| Nl] > = W{W'‘} 



(2.43) 


(2.44) 


where iV* are the one-dimensional Lagrangian quaxlratic shape functions and {W^} and {Wy} 
axe the vectors of the nodal values of the weight functions for a typical boundary element 
shown in Fig. 2.3. Further and ty are approximated as 


iV| Nl] 


( tl 


bx 

2 

bx 

^bx 


= WK4} 


(2.45) 
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and 



’iVf Nl]< 


^ J 

1 


= mH) 


where {t^} and {ty} are the vectors of the nodal values of the traction. 


(2.46) 


2.2.6 Finite Element Equations 


Using the approximations (2.37) and (2.38), the strain rate components and the corre- 
sponding derivatives of the weight functions can be expressed in the vector forms as given 
below. 

The strain rate vector {e} and the vector {e(W)} become 




{£} = [B]{V‘} 


(2.47) 

and 


{«(W")} = 

iB]{wn 


(2.48) 

where 

1 

0 

1 

dx 

0 ‘ 



1 

[5] = 

0 ^ 

^ &y 

0 

a^Q 

(2.49) 



1 ^ 1 ^ 
V2 dy 7 ^ ax 

1 dNa 
U 2 dy 

1 

72 dx 



By substituting the above vector expressions along with the boundary approximations 
(eqn. 2.43 to 2.46) the integral (eqn. 2.33) becomes 

(2-50) 


where 


= ( {SI 1 

= { {SI ) 

IK^] = - l^^[N,f{mflB]dx<ly 


(2.51) 

(2.52) 

(2.53) 

(2.54) 
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[^Ip] - ~ I, [^f{'^}Wp] dx dy 

J 


(2.55) 


[K"J = 2{j.[Bf[B] dx dy + l^j (v^[N,f^ + Vy[N,f^^ dx dy (2.56) 

{fx} = lj^bfW[i]dr (2.57) 

{fy} = {Nt}[i] dT (2.58) 

Here, NEM is the number of area elements, NBEX is the number of boundary elements 
on Fi and NEBY is the number of boundary elements on r 2 . A® denotes the domain of a 
typical area element while Fj and F^ are the domains of typical boundary elements on Fi 
and F 2 , respectively. 

For the purpose of numerical evaluation, the variables of the area integrals (ij.s^o 2 ^rt) 
are transformed to the natural coordinates ( 0 , 77 ) using the following transformation: 

[ {...)dxdy=[ f {...)\J\dC dr} (2.59) 

J J-l 

where \J\ is the determinant of the elemental Jacobian matrix and is given by 

dx dx 
dC dp 

|J| = (2.60) 

^ in 
L 5C 577 J 

Here, the co-ordinates (x,y) are approximated using the biquadratic shape functions. 

Similarly the boundary integrals are transformed to integrals over the interval (-1,1) by 
the relation 

l^{...)dr = j*' {...)\Jt\d<: (2.61) 

where \Jb\ is the Jacobian for the boundary element and is given by 

(2.62) 

Along the boundary, the co-ordinates (x, y) are approximated using 1 -D quadratic shape 
functions. All the elemental matrices are evaluated using the 3X3 Gauss quadrature. Sim- 
ilarly the elemental vectors (2.57 and 2.58) are evaluated using 3 point Gauss quadrature. 

The assembly of elemental coefficient matrices into the global coefficient matrix is done 
by transforming the elements corresponding to a local degree of freedom to positions of 
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corresponding global degrees of freedom in the global coefficient matrix. Similar procedure 
is followed for the assembly of global boundary matrix. The assembled finite element 
equation can be written as 

{Wf[K\{X} = {irf {f } (2.63) 

where {W} is the global vector of nodal values of weight functions, [K\ is the global 
coefficient matrix, {X} is the global vector of nodal values of pressure and velocities and 
{F} is the global right side vector. Since {W} is arbitrary, the final expression is given by 

[K]{X} = {F} (2.64) 

The above equation has to be solved after substitution of the boundary conditions. The 
procedure for the application of boundary conditions is discussed in the next section. 


2.2.7 Application of Boundary Conditions 

The natural boundary conditions on the tool-work interface can be applied without eval- 
uating the boundary matrices {f^} and {fy} in the manner described below. 

The Coulomb friction condition (2.23)at a typical node on the tool work interface (Fig. 
2.3) is given by 

\ts\ = m (2.65) 

After expressing tg and tn in terms of and ty at that node and rearranging the terms, 
the above equation becomes 

tx{cos<f)- f kg sin (j)) - ty {sin (j> + f kg cos<j)) = 0 (2.66) 


where kg is -1 before neutral point and +1 after neutral point. In the present work, (f) is 
assumed to be constant over a boundary element and equal to the angle at the middle node 
of the element because it does not change much over a element. Then the above expression 
for an element can be written as 


{1 - f kg tm(l>) - 


(tan (p + fks) = 0 


(2.67) 


After multiplying above equation by 

/(‘{iV»y{JV»}|4|dC ' (2-68) 
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it becomes 


(2.69) 


{/x}(l - I kstaja<j>)- {/;}(tan0 + / A:,) = {0} 

The above equation holds good at all the nodes of the element. At the middle node 
say ’k’ (global node number), there is no contribution from the neighboring elements and 
therefore, in terms of global right side vector, above equation can be written as 

+ / fc.) = 0 (2.70) 

where dp is the total number of pressure variables. At the end node also, the above eqn. 
(2.70) holds good provided contribution of the adjacent element is added to it. Thus, the 
procedure for applying the condition (2.65) at a boundary node ’k’ is as follows: 

• Replace {dp + 2k — l)th row of global coefficient matrix [K] by the following con- 
dition: 

{l-fkstaxi(l)) times {dp + 2k- l)th row of [AT] -(tan ^6+//:,) times {dp + 2k)th row of [K]. 

• Make the (dp -f 2A: - 1)‘^ row of global right hand vector {F} as zero. 

The velocity boundary condition at the node is applied replacing {dp -\- 2k)th row of [K] 
by the eqn. (1.22). 

• Make {dp -I- 2k)th row of [K] zero. 

• Set {dp + 2k,dp + 2k — l)th element of [K] to tan ^ and {dp + 2k, dp + 2k)th ele- 
ment of [K] to 1. 

• Make {dp + 2k)th row of {F} to zero. • 

2.2.8 Estimation of Neutral Point 

In order to apply the natural boundary condition on the roll-strip interface, the location 
of the neutral point has to be known. In this formulation, the neutral point is found by 
minimizing total power with respect to the position of the neutral point. 

A generalized upper bound theorem [Collins, 1969] is expressed by, 

^ tiV*ds < l^a*je^jdy - tiV*ds (2.71) 

where Sp is the part of surface where some or all of the traction components are specified 
and is the remainder of) the surface. Here, is the boundary CD while Sp includes 
all other boundaries. It is assumed that on parts AB and EF, instead of the velocity 
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component Vx^ it is the stress component tx which is specified with values equal to back 
and front (tf) tensions respectively. Further, is any piece-wise continuous velocity 
field defined over volume V with corresponding strain-rate field and a*j is the plastic 
stress field related to e*j through the flow rule. Here u* is a volume preserving velocity 
field, but does not have to be a kinematically admissible field. But in our formulation, v* 
is so chosen that u* is satisfies all the velocity boundary conditions also. Further, v* Cctn 
be expressed as the sum of roll velocity (Vr) and relative sliding velocity of strip (u; - Vr). 
Therefore, the above inequality becomes, 

tsVRds < a:^€*jdV + ^ tbv:ds - tfV*Js + - Vnjds. (2.72) 

where u* and v* are the normal and tangential components of v* on interface. The integral 
on the left side of the inequality is the actual rolling power and the integrals on the 
right side are plastic power, power due to tensions and Motion power respectively, being 
based on the assumed velocity field. The computations of plastic power and power due to 
tensions pose no problem, whilest for the computation of friction power (last integral of the 
inequality), actual shear traction distribution is required which in not known beforehand. 
Avitzur [1964] has assumed its magnitude to be /(o’y)o/V3) where / is a Motion factor. 
In the present model, is obtained from eqn. (2.23) in which tn is obtained from the 
solution of continuity and momentum equations, at assumed location of neutral point. 
This tn is found to be approximately close to the actual except near assumed neutral 
point. However, since the relative velocity between the roll and strip is very small near 
the assumed neutral point, the error in the evaluation of friction power is negligible. Thus, 
all the three components of power are computed for different assumed positions of neutral 
point and the position which minimizes the sum of these is treated as the correct position 
of neutral point [Dixit and Dixit 1996]. 

2.2.9 Solution Procedure 

After imposing the boundary conditions as discussed in the above section, the final matrix 
is solved iteratively by Householder method [Bathe, 1990]. The Householder method is 
used because the resulting matrix of the mixed formulation is ill-conditioned. For the first 
iteration, the value of ^ is obtained from the upper bound velocity field. For the inlet 
and outlet regions a suitable guess is used. This value is updated by using the velocity 
field of previous iteration. For each assumed position of neutral point, the iterations are 
continued until the nodal velocity and pressure values converge within 0.1% between the 
successive iterations for finding out corresponding power. The position of neutral point 
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corresponding to the minimum power is considered as the correct one. In the process of 
numerical simulation, the value of // becomes very large outside the deformation zone as 
the generalized strain rate is near to zero in those regions. This high value of ^ makes the 
coefficient matrix more ill-conditioned. This difficulty is overcome by prescribing a suitable 
cut off value for (i. 


2.3 Formulation for Strain Hairdening 


To updated the value of /i, in the rigid-plastic formulation, the generalized stress stress a 
need to be updated. In general is a function of e for the metals which are cold rolled. 
Therefore, the procedure to obtain the distribution of e over the domain is presented in 
this section. 

The generalized strain e is obtained by integrating the generalized strain rate I along 
path line (eqn, 2.6) and is given by 


= fut 

Jo 


(2.73) 


Alternately, one can convert the above equation into a differential equation [Avitzur and 
Talbert, 1989] as 

f (2.74) 


de 

dt 


where ^ is the material time derivative. For plane-strain steady-state problem, the above 
equation can be written as 

(2.75) 


de 


de 


Since the deformation at the beginning of the inlet zone negligible, the initial condition for 
this differential equation is 

€ = 0 on AB (2.76) 

Galerkin formulation of eqn. (2.75) over an element gives the following integral: 

de 


^NEM Vy^ - f)l^f dxdy^O (2.77) 

where Wf is the weighting function. The above equation can be written in matrix form as 


^NEM 


r 

" 

ff 1 


J Ae 

{V. V,} < 

1 dz 1 

^ - e 

J /I 


1 ^ J 

- 


Wg dx dy = 0 


(2.78) 
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In the present work, the generalized strain and the corresponding weight function are 
approximated by the same shape functions as used for the velocity approximation. There- 
fore, 

e = {N}{^} ( 2 . 79 ) 

W, = {N}{Wl} (2.80) 

Here vectors {?} and {W^} contain the elemental nodal values of e and Wj. The vector 
containing the derivatives in eqn. (2.78) above can be written as 


di 

'Sx 

. dy 


= PKi*} 


(2.81) 


where 



r Mi 

Mz 

Ma. 1 


dx 

dx 

dx 

[H] = 

dNi 

Mz 



- dy 

dy 

dy J 


(2.82) 


After substituting the approximations for generalized strain, the weighting function, 
the vector containing the strain derivatives and velocity vector, the final expression can be 
written as 


X^NEM 

^e=l 


{wiflKtm 


_ y^NEM 




(2.83) 


where 

[If/] = / [iv.l’'p]) dx dy (2.84) 

Jn 

{fn = I ({Nfi) 'll <i!/ (2-85) 

Jn 

Here, [ATj] and {//} are the elemental coefficient matrix and the right hand side vector 
respectively for the evaluation of generalized strain. Since is arbitrary, the final 

global matrices can be written as 


lii:,l{£-} = {F,} (2.86) 

where [K{\ is the global coefficient matrix, {e} is the global vector of the nodal values of e 
and {Fj} is the global right side vector. 

The generalized strain distribution obtained by solving the above equation is used to 
update the value of generalized stress a. This a is used to update the value of n present 
in the matrices of the rigid-plastic finite element formulations. 
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2.4 Damage Formulation 

2.4.1 Formulation for Critical Damage Criterion [Dhar et al., 
1996] 


In the present work, one of the criteria used for prediction of central burst is the critical 
value of the damage variable. For evaluation of the damage, the continuum damage me- 
chanics model employed by Dhar et al. [1996] is used. The model has been described in 
section 1.2. of chapter 1. For the sake of completeness, some salient features of the model 
have been reproduced below. 


Damage represents surface discontinuities in the form of micro-crack or volume discon- 
tinuities in the form of micro- voids. They assumed that the discontinuities are spread in an 
isotropic way. Then the damage variable is considered as a scalar. In the model proposed 
by Dhar et al. [1996], D is defined as the area void fraction at a point. 

The material behavior of a damaged material is influenced by the damage i.e., the con- 
stitutive equation of such a material depends on the damage variable. For some materials 
like high carbon steel, the damage level at failure are quite low, of the order of 0.05 [Dhar 
et al., 1996]. Therefore, for such materials the damage variable does not affect significantly 
either the deformation field or the stress field. As a result, for such materials, it is possible 
to decouple the damage evaluation law from the constitutive relation. In that case, one 
can first determine the deformation and stress fields using the constitutive relation of the 
undamaged material and then evaluate the damage with the help of damage evaluation 
law. 


Dhar et al. [1996] proposed a damage evaluation law based on the damage growth 
model of Lemaitre [1985] and the experimental results of Ley Roy et al. [1981] and is given 
by 

^ = Cd^+ (oi + a2D){-Y)k (2.87) 

at 

where -Y, the variable conjugate to D, is given by 




( 2 . 88 ) 


+ 3(1 - 2v){^) (2.89) 

The model incorporates the damage growth due to both the void nucleation as well as the 
void growth. 
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For plane-strain problem, the above equation can be written as 


,,dD ^ . 

- {ai + a2D){-Y)e = 0 


(2.90) 


Since the danaage growth at the beginning of the inlet zone is insignificant, the initial 
condition for this difi'erential equation is 

D = 0 on AB (Fig. 2.1) (2.91) 

Galerkin formulation of the eqn. (2.90) over an element gives the following integral: 

- (fli + a 2 D){-Y)e\ Wd dxdy = 0 (2.92) 

where is weighting function. The above integral can be written in the matrix form as 


^NEM f 

Ja^ 


{V. V,} 


I m \ 

dx 


dJD 


\-a2D[-Y)k-{CDA-a,{-Y))k 


Wodxdy = 0 (2.93) 


In the present work, the damage variable D and the corresponding weight function Wd 
are approximated by the same shape functions as used for the velocity approximation. 
Therefore, 

D = (2.94) 

Wd = {N}{WI,} (2.95) 

Here {D^} and {W^} contain the nodal values of D and Wd- The vector containing the 
derivatives of D can be written as 


dx 

9y 


= [H]{D^} 


(2.96) 


where [H] has already been defined in eqn. (2.82). After substituting the approximations 
for D, Wd, the vector containing the derivatives of D and the velocity vector, the final 
expression can be written as 




(2.97) 


where 


[^d] = - G 2 (-r){iV}^{N}?] dx dy (2.98) 

{/^}= [ {N}'^{CD + ai{-Y))kdxdy (2.99) 

J 



Here [Kj^] and {/j^} are the elemental coefficient matrix and right side vector respectively 
for the evaluation of damage D. Since is arbitrary, the final global matrices can be 

written as 

[Kd]{D} = {Fd} ( 2 . 100 ) 

where [Kd] global coefficient matrix, {-D} is the global vector of the nodal values of 
D and {Fd} is the global right side vector. 

While solving the above equation, Y is calculated from the deviatoric stress field and 
the pressure field of the rigid-plastic solution of eqn. (2.64). 


2.4.2 Formulation for Oyane’s Criterion [Oyane et al., 1980] 


In this section, Galerkin finite element formulation for Oyane’s fracture criterion is pre- 
sented for prediction of defects in rolling, which is proposed by Oyane et al. [1980]. The 
model has been described in section 1.2. of chapter 1. 

Oyane et al. [1980] used the plasticity theory for porous materials to propose a duc- 
tile fracture criterion indicating that micro crack initiates whenever the volumetric strain 
reaches a material dependent critical value C. They proposed a empirical criterion which 
is given by 

where e/ is the generalized strain at which fracture occurs, C and A are material constants. 
This integral can also be written as 


and also. 




dC , ^m\L. 

= + 


For plane-strain problem, the above equation can be written as 

T/ T/ Al I 


( 2 . 102 ) 

(2.103) 

(2.104) 


Since the Damage growth at the beginning of the inlet zone is insignificant, the initial 
condition for differential equation is 

C = 0 on AB' (2.105) 

Galerkin formulation of the eqn. (2. 104) over an element gives the following integral: 
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where Wc is weighting function, the above integral can be written in the matrix form 3 S 


^NEM 


L 


{K y,} 


dc 

dx 

dc 

dy 



Wc dxdy = 0 


(2.107) 


In the present work, the damage variable C and the corresponding weight function Wc 
are approximated by the same shape functions as used for the velocity approximation. 
Therefore, 

= (2.108) 

Wc = {N}{W^} (2.109) 

Here {C**} and {W^} contain the nodal values of C and Wc- The vector containing the 
derivatives of C can be written as 


( ^ \ 
dx 

< > 
< dy - 




( 2 . 110 ) 


where [H] has already been defined in eqn. (2.82). After substituting the approximations 
for C and Wc, the vector containing the derivatives of C and the velocity vector, the final 
expression can be written as 


mfyK’cMc’} = ( 2 . 111 ) 

where 

lK‘c] = dx dy (2.112) 

Here [K^] and {fc} are the elemental coefficient matrix and right side vector respectively 
for the evaluation of damage C. Since {W^}'^ is arbitrary, the final global matrices can be 
written as 

[Kc]{C} = (2.114) 

where [Kc] is the global coefficient matrix, {C} is the global vector of the nodal values of 
C and {Fc} is the global right side vector. 

In order to apply the Oyane’s criterion, the availability of the material dependent 
constant A and C is necessary. With this aim Alberti et al. [1993] have been carried 
out some upsetting tests on cylindrical specimens of aluminium alloy (UNI 3571). They 
obtained values of A and C, they are 0.24 and 0.348 respectively. 
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2.5 Secondary Quantities 


Prom the above formulations one can obtain the primary variables only, namely the nodal 
velocities, pressure, deviatoric stresses etc. The procedure for obtaining the secondary 
quantities from them is described in this section. 


Roll Torque (T) 


The roll torque is calculated from the relationship 

PR 


T = 




(2.115) 


where the total power P consists of the following three parts. 


Power required for plastic deformation (Pp) 

The power dissipated due to plastic deformation is given by 

■Pp = / c^ij^ij dx dy (2.116) 

J A 

Power required to overcome friction at the roll-strip interface (Pj) 

The power dissipated due to friction is given by 

Pf= tMAVs\ds (2.117) 

^0 

where is the relative velocity with respect to the velocity of neutral point along the 
interface and I is the arc of contact. The shear stress t, is calculated from eqn. (2.23). 


Power due to tensions (Pj) 

The power required in the presence of front and back tensions is 

Roll pressure of interfacicd normal stress (tn) 

At a typical boundary node shown in Fig. (3), the normal pressure can be expressed as 

tn = tx rix + ty riy (2.119) 
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where riy — cos(^ and — sin(^. now by using expression ti = cTy-nj the above equation 
can be written as 

in = (^xx^x^ + (yyyTiy^ + 2axynxny ( 2 . 120 ) 

To get the stress value at the interface, first the stresses are calculated at 3X3 Gauss 

points in the element and then extrapolated to 3 points on the interface. 

§ 

Roll force (Fr) 

The roll force is the component of the resultant of interfacial stresses which is given by 




Figure 2.1: The Domain and the boundary conditions (velocity and traction) for plane- 
strain rolling. 



Figure 2.2: Typical area element. 



C 

Figure 2.3: Typical boundary element. 
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Chapter 3 

Results and Discussions 


3.1 Introduction 

The finite element model presented in chapter 2 is implemented using Fortran language and 
applied to analyze plane-strain rolling process. The results are obtained using three metals. 
The metals considered being aluminium, copper and steel. These metals are chosen, so the 
comparison with experimental work is facilitated. 

In section 3.2, fracture criteria namely hydrostatic stress criterion [Reddy et al., 1996], 
critical damage criterion [Dhar et al., 1996] and Oyane’s criterion [Oyane et al., 1980] are 
presented. 

To test the predictions (pressure, force and torque) of the present model, its predic- 
tions are compared with the experimental results of Al-Salehi et al. [1973] and Shida and 
Awazuhara [1973] and also with the finite element results of . Dixit [1997] and are 

presented in sections 3.3.1 and 3.3.2.' In section 3.3.3, the predictions of HSC are compared 
with the available experimental results of Turckzyn and Malinowski [1996] for alhgatoriug. 

Comprehensive parametric study is carried out to study the influence of process vari- 
ables on the hydrostatic stress and damage distributions [Dhar et al., 1996; Oyane et al., 
1980] and are presented in sections 3.4.1 to 3.4.3. Parametric study for studying the effect 
of process variables on the boundary between safe and unsafe zone for fracture is carried out 
using fracture criteria namely hydrostatic stress criterion [Reddy et al., 1996] and critical 
damage criterion [Dhar et al., 1996] and the same is presented in section 3.4.4. 
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3.2 Fracture Criterion 


DifFersnt fracture criteria used to predict the fracture initiatioii are presented in this section. 
Three criteria namely hydrostatic stress criterion [Reddy et al., 1996], critical damage 
criterion [Dhar et al., 1996] and Oyane’s criterion [Oyane et al., 1980] are used to predict 
the central burst and only HSC is used to predict the split-ends in rolling. 

3.2.1 Hydrostatic Stress Criterion (HSC) [Reddy et al., 1996] 

If the process conditions are such that hydrostatic stress is tensile in the deformation zone, 
then micro-voids first get nucleated, then they grow and finally coalesce with one another 
leading to micro-crack initiation. If one can maintain the compressive hydrostatic stress in 
the deformation zone, then the micro-voids either do not nucleate or the existing micro- 
voids remain closed and, thus, there is no possibility of micro-crack initiation. The work 
of Clift et al. [1990] shows that the value of hydrostatic stress at fracture initiation in 
extrusion is close to zero. Based on these observations Reddy et al. [1996] proposed the 
hydrostatic stress criterion to predict the central burst in extrusion, and the same is used 
in the present work to predict the central burst in rolling. 

Criterion for Central burst 

Whenever hydrostatic stress at a point on the centerline in the deformation zone becomes 
zero and is compressive elsewhere, there is fracture initiation leading to central burst. 

Criterion for Alligatoring 

When the damage is severe it leads to alligatoring or split-end defect. Alligatoring also 
nucleates as a micro- void, leading to micro-crack, forming along the central plane of the 
deformed material (the plane connecting the axis of symmetry of the rolls) . As the rolling 
proceeds two ’’halves” of the material separates from each other and alligatoring occurs. It 
will occur when the tensile hydrostatic stress at the axis of symmetry reaches the order of 
yield stress of the material. Based of this hydrostatic stress criterion to predict split-end 
is proposed as: 

’’Whenever the tensile hydrostatic stress value at a point on the centerline in the de- 
formation zone reaches the order of yield stress, the split-end defect will occur . 
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3.2.2 Critical Damage Criterion [Dhar et al., 1996] 

Dhar et al. [1995] reported that for SAE 1090 steel, micro-crack initiation arises whenever 
the value of damage variable D reaches the critical value (D^) of 0.05. The same is used 
in present work to predict the central burst initiation in rolling. 

3.2.3 Oyane’s Criterion [Oyane et al., 1980] 

Alberti et al. [1993] reported that for aluminium alloy (UNI 3571), fracture occurs when- 
ever the value of damage variable C in Oyane’s criterion [Oyane et al. 1980] reaches the 
critical value (Cc) of 0.348. 

3.3 Validation 

3.3.1 Comparison of Pressure Predictions 

Pressure predictions of the present model are compared with the experimental results 
of Al-Salehi et al. [1973]. Figures (3.1) and (3.2) shows the variation of roll pressure 
with the angle of contact for copper and aluminium respectively. It can be seen that the 
predicted results are in reasonable agreement with the experimental results of Al-Salehi et 
al. [1973]. The general shapes of the experimental and predicted distributions are similar, 
there being differences in magnitude. The difference in magnitude can be attributed to the 
roll deformation which is neglected in the present work. 

Pressure predictions of the present model are compared with the finite element results 
of Dixit , [1994] (Figures (3.1) and (3.2)). It can be seen from the figures (3.1) and 

(3.2) that predictions of both the models are in good agreement. However shift in neutral 
point towards the exit has been observed in case of aluminium. The discrepancy between 
the two results can be attributed to roll deformation and variation of friction coefficient, 
which is considered by . , Dixit [1994]. 

3.3.2 Comparison of Force and Torque Predictions 

Force and torque predictions of the present model are compared wnth the experimental 
results of Al-Salehi et al. [1973] for aluminium and copper and Shida and Awazuhara 
[1973] for steel. Figures (3.3) to (3.5) shows the variation of roll separating force with the 
reduction in area for aluminium, copper and steel respectively. Figures (3.6) to (3.8) shows 
the variation of roll torque with the reduction for aluminium, copper and steel respectively. 
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It can be seen from the figures 3.3 to 3.8, at the lower reductions the discrepancy between 
the predicted and experimental results is low. At higher reductions the discrepancy between 
experimental and predicted results is high. This can be attributed to increase in roll 
deformation with reduction. In case of aluminium (Figures 3.3 and 3.6) and copper (Figures 

3.4 and 3.7) the discrepancy at higher reductions also very less compared to steel (Figures 

3.5 and 3.8), this is mainly due to the roll deformation which will be higher in case of steel. 
In the present model roll deformation is not considered. 

Roll force and torque predicted in the present model are also compared with the finite 
element predictions of Dbcit [1997] (Figures (3.3) to (3.8)). There is a good 

agreement between predicted and finite element results of Dixit [1997]. However 

there is a discrepancy in magnitude this can be attributed to the roll deformation and 
variation in friction which is not considered in the present work. 

3.3.3 Validation of HSC 

Reddy et al. [1996] proposed hydrostatic stress criterion to predict central burst in extru- 
sion. They compared the prediction of HSC with the available experimental results and 
showed that there is a good agreement between HSC predictions and experimental results. 
The same HSC is used in the present work to predict central burst in rolling. Critical 
literature survey carried out to find the experimental results to compare the prediction of 
HSC for central burst but, the attempt was unsuccessful. 

To establish the validity of the HSC predictions for spUt-ends/ alligatoring of the present 
work, these results are compared with the available experimental results of Turckzyn and 
Malinowski [1996]. They carried out experiments using an aluminium alloy (6061-T6) of 
rectangular bars of 6.35, 7.94, 9.35 and 12.7 mm of height with 25.4 mm of width. They 
used two pairs of rolls having radii 21 and 67 mm, to obtain results at eight values (0.095, 
0.118, 0.139, 0.189, 0.302, 0.378, 0.44 and 0.605) of parameter (ho/i?o). These experimental 
condition may not be satisfying plane-strain deformation conditions but, they compared 
these results with their plane-strain upper bound results. The material properties for 
annealed aluminium alloy as give by Turckzyn and Malinowski [1995] are 

Generalized stress (d) — 24.83[1 -f 57.74 e] MPa (3.1) 

Generalized stress at zero plastic strain (ao) = 43 MPa (3-2) 

Avitzur et al. [1988] also used the same parameter (/lo/ Ro) to obtain the safe and unsafe 
zone boundary for central burst as well as split-ends. Hence, from these works [Avitzur et 
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al.) 1988, Turckzyn and Malinowski, 1996], it can be concluded that safe and unsafe zone 
for defect is independent of {ho) and (-Ro)) aiid depends on the parameter {ho/Ro). 

Turckzyn and Malinowski [1996] uSed roll of 67 mm radius, to obtain results at first 
four values of parameter {ho/Ro) { 0.095, 0.118, 0.139 and 0.189). It can be seen from 
the Fig. 3.9 that they carried out experiments in the range of 0.03 to 0.16 reductions. 
Except few stray points, defect is not occurred in this range. For remaining four values of 
parameter {ho/Ro) (0.302, 0.378, 0.44 and 0.605) using roll radius of 21 mm, they carried 
out experiments in the range of 0.05 to 0.27 reductions. Fig. 3.9 shows that the defect 
occurred upto 0.16, 0.18 and 0.24 reductions at ho/Ro values of 0.302, 0.378 and 0.44 
respectively. 

In the present Avork, to establish the validity of the hydrostatic stress criterion, its pre- 
dictions for split-ends are compared with experimental results of Turckzyn and Malinowski 
[1996]. It is assumed that alligatoring occurs whenever the ratio of hydrostatic stress to 
material yield stress ratio {cm/oo) is within the limits of -0.85 and -1.15. Occurrence of 
split-ends is predicted at values of ho/Ro from 0.095, 0.118, 0.139, 0.189, 0.302, 0.378 and 
0.44 at the same conditions as used by Turckzyn and Malinowski [1996]. In the present 
work friction coefficient is taken as 0.1. To predict the split-end occurrence at values of 
ho/Ro from 0.095 to 0.189, using roll radius 67 mm, executions are carried out. It is ob- 
served that from reduction 0.03 to 0.17 there is no occurrence of defect. Using roll radius 
21 mm executions are carried out at three values (0.302, 0.378 and 0.44) of parameter 
ho/Ro and defects are not predicted above the reductions 0.19, 0.25 and 0.29 respectively. 
Although, the thickness to width ratio used for carrying out experiments by Turckzyn and 
Malinowski [1996] is not in the plane-strain range, they are compared with the predictions 
of the present work which assumes the deformation to be plane strain. From the above ob- 
servations, it can be seen that the predictions of the present work using HSC for split-ends 
are in good agreement with the experimental results of Turckzyn and Malinowski [1996]. 

It is observed that keeping {ho/Ro) parameter value constant, if the values of ho and Ro 
are changed, the hydrostatic stress distribution remains same but the value of compressive 
hydrostatic stress is varying. This can be seen from Figures 3.10 (a) and (b). Hence, 
the boundary between safe and unsafe zone is not independent of ho and Ro and also not 

dependent on {ho/Ro)- 
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3.4 Parajmetric Study 


Parametric study is carried out to study the influence of process variables (reduction in 
area, friction and roll radius) on the hydrostatic stress distribution, damage distribution of 
critical damage criterion [Dhar et al., 1996] and damage distribution of Oyane’s criterion 
[Oyane et al., 1980]. Parametric study for studying the effect of process variables on the 
safe and unsafe zone for fracture is carried out using fracture criteria (hydrostatic stress 
criterion [Reddy et al., 1996] and critical damage criterion [Dhar et al., 1995) same 

is presented in section 3.4.4. 

3.4.1 Hydrostatic Stress Distribution 

The results in this section are obtained by choosing SAE 1090 steel, the material properties 
for SAE 1090 steel as given by Le Roy et al. [1981] are 


Generalized stress (a) 

= 1115 (e)°-^® MPa 

(3.3) 

Generalized stress at zero plastic strain (co) 

= 464 MPa 

(3.4) 

Modulus of Elasticity (E) 

= 210 GPa 

(3.5) 

Poisons ratio (v) 

= 0.3 

(3.6) 


Hydrostatic stress plays a major role in ductile fracture. Therefore, for prediction of ductile 
fracture, it is necessary to know the effect of process variables on the hydrostatic stress 
distribution. Figure 3.11 shows the hydrostatic stress distribution for a typical set of 
process conditions (Material: SAE 1090 steel, % r = 10, f = 0.075, R = 65 mm). It is 
observed that the hydrostatic stress is compressive everywhere in the deformation zone, but 
becomes tensile at the boundary in a small region near the exit of the roll-work interface. 
Figure 3.12 shows that the hydrostatic stress distribution remains same with increase in 
reduction in area but the value of compressive hydrostatic stress is increasing with increase 
in reduction. When the friction is decreased (Fig. 3.13), the distribution remains same 
with an increase in compressive stress value in the deformation zone. With increase in roll 
radius (Fig. 3.14) the compressive hydrostatic stress is decreasing. This can be attributed 
to increase in homogeneity of deformation with roll radius. 
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3.4.2 Damage Distribution [Dhar et al., 1996] 

The results in this section are obtained by choosing SAE 1090 steel as the work piece 
material. This material is chosen because of the availability of data that axe required for 
applying the fracture criterion of Dhar et al. [1996]. The material properties of SAE 1090 
steel are given in above section. 

Figure 3.15 shows the damage distribution for a set of process parameters (Material: 
SAE 1090 steel, % r = 10, f = 0.075, R = 65 mm). It is observed that value of the damage 
variable is increasing along the plane of symmetry and reaching maximum at exit of the 
roll- work interface. Fig. 3.16 shows that there is no change in damage distribution with 
increase in reduction but the value of damage increasing in the deformation zone. This 
can be attributed to the increase in inhomogeneity of the deformation with increase in 
reduction. Fig. 3.17 shows that there is no significant change in damage distribution with 
decrease in friction coefficient. However, the maximum value of damage variable is slightly 
decreasing with decrease in friction coefficient. Fig. 3.18 shows the variation of damage 
distribution with roll radius. With the increase in roll radius the max value of damage in the 
deformation zone is decreased, this is mainly due to increase in homogeneity of deformation 
with roll radius. In all these (Figures 3.15 to 3.18) damage distribution contours damage 
variable is multiplied by 10. 

3.4.3 Damage Distribution [Oycine et al., 1980] 

The results of this section axe obtained by choosing aluminium alloy (UNI 3571) as the 
workpiece material. This material is chosen because of the availability of data that axe 
required for applying the Oyane’s criterion [Oyane et al., 1980]. Properties of this alloy 
given by [Alberti et al., 1993] 

Generalized stress (a) = 574 (e) MPa (3.7) 

Generalized stress at zero plastic strain (cto) = 300 MPa (3.8) 

Figure 3.19 shows the damage distribution for a typical set of process parameters (% 
r = 10, f = 0.075, R = 65 mm). It is observed that value of the damage variable is 
increasing along the plane of symmetry and reaching maximum at exit of the roll-work 
interface. It is observed that with the increase in reduction the ma xim u m damage value 
also increasing along the plane of symmetry (Fig. 3.20). This can be attributed to the 
increase inhomogeneity of deformation. The maximum value of damage value is decreasing 
with decrease in friction (Fig 3.21), this can be attributed to the decrease in inhomogeneity 
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of deformation. With increase in roll radius the maxiinuin value of the damage is decreasing 
(Fig 3.22). This can be attributed to the increase in deformation length with roll radius 
there by decrease in inhomogeneity of deformation. 

3.4.4 Safe and Unsafe Zones for Defects 

3.4.4. 1 Hydrostatic Stress Criterion 

Predictions of this section are obtained by using SAE 1090 steel, the properties of this 
material are given in section (3.2.3). 

Alligatoring 

Fig. 3.23 shows the boundary between safe and unsafe zone for perfectly-plastic and strain 
hardening material. It can be observed that for strain hardening material the safe zone 
is less compared to perfectly plastic material. This trend is similar to that of Turckzyn 
and Malinowski [1996]. Fig. 3.24 shows the variation of safe and unsafe zone with respect 
to variation in coefficient of friction. As the coefficient of friction increases the boundary 
between safe and unsafe zone is moving towards higher reduction, i.e., reduction in safe 
zone. This can be attributed to decrease in hydrostatic stress with increase in friction. 
Where as the zone is moving towards less reduction with increase in roll radius (Fig. 3.25), 
this can be attributed to the increase in deformation length with increase in roll radius 
there by increase in homogeneity of deformation. A typical hydrostatic stress distribution 
at split-end condition is shown in Fig. 3.26. 

Central Burst 

Figures 3.27 and 3.28 shows the change in boundary between safe and unsafe zone with 
change in roll radius and coefficient of friction respectively. It can be observed that safe 
zone is increased with decrease in friction or increase in roll radius. These effects are 
similar to that of alligatoring observed in above section. Fig. and 3.29 shows the typical 
hydrostatic stress at central bursting condition. 

3. 4. 4. 2 Critical Damage Criterion [Dhatr et al., 1996] 

Predictions of this section are obtained by using SAE 1090 steel, the properties of this 
material are given in section 3.4.2. 


' ft 
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Figures 3.30 and 3.31 shows the changes in safe and unsafe zone boundaries with change 
in roll radius and friction coefficient respectively. Figures 3.30 shows that with increase in 
roll radius the boundary between safe and unsafe zone is moving towards lower reductions. 
Figures 3.31 shows that with friction coefficient the boundary between safe and unsafe zone 
is moving towards higher reduction. These predictions shows the similar trend as the HSC 
predictions shown in above section. A typical damage distribution at fracture conditions 
is shown in Fig. 3.32. In Fig. 3.32 damage variable is multiplied by 10. 

It can be observed that predictions of HSC for central burst is always above the criticed 
damage criterion. Hence it can be said HSC is more conservative as mentioned by Reddy 
et al. [1996]. 
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Figure 3.1: Comparison with experimental and finite element results for roll pressure [Ma- 
terial: aluminium, a = 50.3 (1 + e/0.050) MPa, R = 79.375 mm, ho = 6.274 mm, f = 
0.1, % r = 14.17] 



Figure 3.2: Comparison with experiment^ and finite element results for roll pressure [Ma- 
terial: copper, a = 70.3 (1 4- e/0.022) MPa, R = 79.375 mm, ho = 6.35 mm, f = 0.06, 
% r = 14.80] 
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Figure 3.3: Roll force comparison with finite element and experimental results [Material: 
aluminium, a = 50.3 (1 + e/0.050) MPa, ho = 6.274 mm, R = 79.375 mm, f = 0.10] 



Figure 3.4: Roll force comparison with finite element and experimental results [Material: 
copper, a = 70.3 (1 + e/0.022) MPa, ho = 6.35 mm, R = 79.375 mm, f = 0.06] 
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Figure 3.5: Roll force comparison with finite element and experimental results [Material; 
steel, a — 324.0 (1 + e/0.052) MPa, Hq = 1 mm, R = 65 mm, f = 0.08] 



Figure 3.6: Roll torque comparison with finite element and experimental results [Material: 
aluminium, a = 50.3 (1 + e/0.050) MPa, ho = 6.274 mm, R = 79.375 mm, f = 0.10] 
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Figure 3.9: Analytical and experimental results of Turckzyn and Malinowski [1996] for 
split-ends 
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Figure 3.10: Hydrostatic stress (^) distribution at {ho/Ro) = 0.1 (a) ho - I mm, fZo - 65 
mm, (b) ho = 2 mm, Rq = 130 mm [Material: SAE 1090 steel, a = 1115 MPa, f = 




Figure 3.11: Hydrostatic stress (^) distribution [Material; SAE 1090 steel, c 
1115 MPa, % r = 10, f = 0.075, R = 65 mm] 



Figure 3.12: Hydrostatic stress (^) distribution [Material; SAE 1090 steel, o 
1115 MPa, % r = 15, f = 0.075, R = 65 mm] 
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Figure 3.13: Hydrostatic stress (^) distribution [Material: SAE 1090 steel, a 
1115 MPa, % r = 10, f = 0.05, R = 65 mm] 



Figure 3.14: Hydrostatic stress (^) distribution [Material: SAE 1090 steel, a 
1115 MPa, % r = 10, f = 0.075, R = 130 mm] 
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Figure 3.15: Damage distribution in plane-strain rolling, critical damage criterion [Dhar 
et al., 1996] (DXIO) [Material: SAE 1090 steel, a = 1115 e MPa, % r = 10, f = 0.075, 
R = 65 mm] 



Figure 3.16: Damage distribution in plane-strain rolling, critical damage criterion [Dhar 
et al., 1996] (DXIO) [Material: SAE 1090 steel, a = 1115 e MPa, % r = 15, f = 0.075, 


R = 65 mm 
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Figure 3.17: Damage distribution in plane-strain rolling, critical damage criterion [Dhar 
et al., 1996] (DXIO) [Material: SAE 1090 steel, a = 1115 e MPa, % r = 10, f = 0.05, 
R = 65 mm] 



Figure 3.18: Damage distribution in plane-strain rolling, critical damage criterion [Dhar 
et al., 1996] (DXIO) [Material: SAE 1090 steel, a = 1115 e MPa, % r = 10, f = 0.075, 
R = 130 mm] 
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Figure 3.19: Damage distribution in plane strain rolling, Oyane’s criterion [Oyane et al., 
1980] (Material: aluminium alloy (UNI 3571), d = 574 MPa, % r = 10, f = 0.075, R 
= 65 mm) 



gure 3.20; Damage distribution in plane strain rolling Oyane's criterion [Oyane etal 
80] (Material: aluminium alloy (UNI 3571), ^ = 574 MPa. % r = 15, f = 0.075, R 




Figure 3.21: Damage distribution in plane strain rolling, Oyaue’s criterion [Oyane et al., 
1980] (Material: aluminium alloy (UNI 3571), o — 574 MPa, % r = 10, f = 0.05, R 
= 65 mm) 



Figure 3.22: Damage distribution in plane strain rolling, Oyane’s criterion [Oyane et al., 
1980] (Material: aluminium alloy (UNI 3571), a = 574 MPa, % r = 10, f = 0.075, R 
= 130 mm) 
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Figure 3.23: The boundary between safe and unsafe zones for perfectly plastic and strain 
hardening using HSC for split-ends [Material: SAE 1090 steel, a = 1115 MPa, R = 
67 mm, f = 0.1] 



Figure 3.24: The boundary between safe and unsafe zones for different friction using HSC 
for split-ends [Material: SAE 1090 steel, a = 1115 e°-^AMPa, R = 65 mm] 
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InitiaJ thickness (ho/Ro) 


Figure 3.25: The boundary between safe and unsafe zones for different roll radius using 
HSC for split-ends [Material: SAE 1090 steel, a = 1115 e MPa, f = 0.08] 



Figure 3.26: A typical hydrostatic stress distribution at split-end condition [Material: SAE 
1090 steel, a = 1115 MPa, R = 65 mm at ho = 3 mm, f = 0.08] 





Figure 3.27: The boundaxy between safe and unsafe zones for different roll radius using 
HSC for central burst [Material: SAE 1090 steel, a = 1115 MPa, f = 0.08] 



Figure 3.28: The boundary between safe and unsafe zones for different friction using HSC 
for central burst [Material: SAE 1090 steel, a = 1115 MPa, R = 65 mm] 
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Figure 3.29; A typical hydrostatic stress distribution at central bursting condition (using 
HSC) [Material: SAE 1090 steel, a = 1115 MPa, /iq = 2 mm, R = 65 mm, f = 0.08] 



Figure 3.30: The boundary between safe and unsafe zones for different roll radius using 
critical damage criterion [Dhax et al., 1996] [Material; SAE 1090 steel, a = 1115 e MPa, 
f = 0.08] 
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Figure 3.31: The boundary between safe and unsafe zones for different friction conditions 
using critical damage criterion [Dhar et al., 1996] [Material: SAE 1090 steel, a = 1115 e 
MPa, R = 65 mm] 



Figure 3.32: A typical hydrostatic stress distribution at central bursting condition (Dhar 
et al., 1996) [Material: SAE 1090 steel, a = 1115 e°-^® MPa, /lo = 4 mm, R = 65 mm, f = 
0.08] ’ 
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Chapter 4 


Conclusions and Scope for Future 
Work 

4.1 Conclusions 

On the basis of results presented in chapter 3, the following conclusions are made. 

• The present predictions of force, torque and pressure axe in good agreement with the 
available experimental [Al-Salehi et al., 1973 and Shida and Awazuhara, 1973] and 
analytical [Dixit and Dixit, 1997] results. 

• The present predictions based on hydrostatic stress criterion for split-ends are in good 
agreement with the available experimental [Turckzyn. and Malinowski, 1996] results. 

• Parametric study on the damage variable of critical damage criterion [Dhar et al., 
1996] and Oyane’s criterion [Oyane et al., 1980] shows that, damage distribution 
remains same but the value of damage variable is increasing with decrease in roll 
radius and increase in friction or reduction. 

• Parametric study on the boundary between safe and unsafe zone is carried out to 
study the influence of process variables. It is observed that the safe zone is increasing 
with increase in roll radius or decrease in friction. The safe zone is higher for non- 
hardening material compared to the hardening material. 

• It is observed that the predictions of HSC are more conservative than the critical 
damage criterion [Dhar et al., 1996]. 
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4.2 Scope for the Future Work 


In the present work, roll deformation, variation of coefficient of Motion along the roll-work 
interface, effect front and back tensions are. not considered and they can be included. To 
validate the fracture predictions made by the present work, experiments can be conducted. 
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